Small Fermi energy, zero point fluctuations and nonadiabaticity in MgB 2 
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Small Fermi energy effects are induced in MgB2 by the low hole doping in the a bands which 
are characterized by a Fermi energy Ep ~ 0.5 eV. We show that, due to the particularly strong 
deformation potential relative to the Ei g phonon mode, lattice fluctuations are reflected in strong 
fluctuations in the electronic band structure. Quantum fluctuations associated to the zero-point 
lattice motion are responsible for an uncertainty of the Fermi energy of the order of the Fermi en- 
ergy itself, leading to the breakdown of the adiabatic principle underlying the Born-Oppenheimer 
approximation in MgB2 even if uip^/E-p ~ O.f — 0.2, where u; p h are the characteristic phonon frequen- 
cies. This amounts to a new nonadiabatic regime, which could be relevant to other unconventional 
superconductors. 
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Four years after the discovery of superconductivity at 
39 K in MgB 2 ,- the superconducting properties of this 
material are still object of investigation. One of the 
anomalous characteristics of MgB2 with respect to con- 
ventional ME superconductors is the remarkable small- 
ness of the Fermi energy Ep associated with the hole- 
like a bands, which are the most involved in the Cooper 
pairing. ME theory holds true indeed only if the Fermi 
energy E-p is much larger than any other energy scale of 
the system. Several analyses indicate however that this 
is not the case of MgB2. LDA calculations, for instance, 
estimate the energy distance between the chemical po- 
tential and the top of the a bands to be ~ 0.4-0.6 eVj^ 
in agreement with ARPES measurements. 4 Estimates of 
the Fermi energy can be inferred also from penetration 
depth measurements giving Ep = 3200 K in MgB 2 ^ and 
even less in Al and C doped compounds^ 

In previous studies we have analyzed some of the effects 
related to a small Fermi energy Ep on different electronic 
and vibrational properties of MgB 2 . One of these effects 
is the breakdown of Migdal's theorem, and consequently 
of the Migdal-Eliashberg diagrammatic theory, which oc- 
curs when the Fermi energy Ep and the phonon energy 
uj p h are comparable (uj p h/Ep < Another one is the 

remarkable anharmonicity of the E 2g phonon mode, re- 
lated to the splitting in energy of the a bands under the 
E 2g distortion^ 

In this paper we want to discuss in greater detail new 
physical consequences of the comparable size of Ep and 
the E 2g splitting energy of the a bands in MgB 2 . In par- 
ticular we review in a critical way the validity of the adi- 
abatic Born-Oppenheimer (BO) approximation for the 
E2 g phonon, which is the most relevant to the supercon- 
ducting pairing. As our main result we show that, due to 
the large quantum fluctuations associated with the zero 
point motion, the BO principle is broken down indepen- 
dently of the ujph/Ep ratio. This novel kind of adiabatic 
breakdown is thus shown to be related to the parameter 
k = gE 2 /Ep, where gE 2g is the electron-phonon matrix 
element which couples the electrons in a bands to the 



E 2g phonon mode in MgB 2 . A simple quantum analysis 
gives k ~ 0.9 — 1, strongly questioning the reliability of 
a Born-Oppenheimer based analysis. 

The Born-Oppenheimer principle is implicitly assumed 
in writing down an effective electron-phonon model start- 
ing from ab-initio techniques. In simple terms, the BO 
principle permits to describe the full many-body wave 
function \l/(r;u) in terms of two separate quantum sys- 
tems: a purely electronic problem, which depends para- 
metrically on the ion variables, and an effective lattice 
problem, where the electronic degrees of freedom have 
been integrated outiiS The total wave function ^f(r; u) 
can thus be written as the product of two "partial" wave 
functions: 

(u)vn(r; [«]), (1) 

where r and u represent the electronic and lattice degrees 
of freedom respectively and (. . . ; [u]) indicates a para- 
metric (not quantum) dependence on the lattice variable. 
The index n identifies the quantum number of the elec- 
tronic state and the label a the phonon cigenstate in the 
effective lattice potential which depends on the electronic 
cloud <p n (r; [u]). 

In the spirit of the BO approximation, the elec- 
tronic ground state wave function ip(r; [u]) is obtained 
as the lowest energy electronic state available for fixed 
u. M. Born and J.R. Oppenheimer showed that Eq. 
||TJ can be considered a good approximation as far as 
(— h 2 /M )V^ip(r; [u]) (M being the atomic mass) is neg- 
ligible. This assumption corresponds to the adiabatic 
hypothesis that the electronic dynamics of (p(r; [u]), is 
much faster than the lattice one. In common metals 
a widely used parameter to evaluate this condition is 
the adiabatic ratio uj p h/Ep, where Ep can be estimated 
in the undistorted u — ground state configuration, 
E-p = Ep[u = 0], and w p h is the lowest energy excita- 
tion of the Xa,n phonon spectrum. However it should be 
noted that, in order for the adiabatic principle to apply, 
the condition (-^ 2 /M)V 2 tp(r; [u]) ~ must be fulfilled 
in all the region of the u-space where x(u) has a sizable 
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weight. This constraint results to be much stronger than 
the conventional adiabatic condition, ajp^/ E-p -C 1, since 
it implies the adiabatic ratio to be negligible for generic 
lattice distortions u, namely lo^/ E-p[u] -C 1, where the 
lattice wave-function |x(u)| 2 provides the range of the 
relevant it-space. As we are going to see, the breakdown 
of this assumption for a relevant range of the Ei g lattice 
distortion u is expected to be responsible for the partial 
failure of the BO approximation in MgB2. 

In order to test the validity of the BO approximation, 
the example of frozen phonon techniques is particularly 
instructive since they are a direct implementation of the 
BO principle. In the frozen-phonon approach the elec- 
tronic structure and the total energy of the system are 
calculated as a function of a static lattice distortion, in 
the spirit of the BO principle, providing ip(r; [u]). The 
frozen phonon lattice potential is therefore used to solve 
the lattice quantum problem and to get the ground state 
lattice wave function x( u )- We shall now apply the same 
technique to MgB 2 , in order to show that, due to the 
small Ep and to the strong deformation potential of 
the E2 g phonon mode, the BO approximation cannot be 
safely applied in this system. 

As mentioned above, small Fermi energy effects in 
MgB2 involve the subset of the almost two-dimensional 
a bands with a corresponding Fermi energy Ep = 0.4- 
0.6 eVj£ defined as the distance in energy between the 
chemical potential and the edge of the a band. Fur- 
thermore, the electron-phonon interaction is essentially 
dominated by the doubly-degenerate Ei g phonon mode 
at q = OfSiii which gives rise to a significant a-o band 
electron-phonon coupling, while the residual a-ir and 7r- 
7r electron-phonon interaction is much weaker^* 1 * The 
effective Hamiltonian of the a bands interacting with the 
E2 g phonon mode at q = can be thus written as£ 



-2y 



2M 



E 2g 
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where u is the lattice displacement of the Ei g phonon 
mode and the sign ± in the Jahn- Teller electron-phonon 
interaction is different for the two a bands. Assuming 
for the moment a perfectly harmonic phonon mode with 
elastic constant a 2 : V(u) = a2U 2 , we can write Eq. (J2J 
in second quantization [u = (h 2 / %M E2 a2) 1 ' i (a + a*)]: 



H = e k 4c k ± g E , 



k 



4 c k(a 



(3) 



where the electron-phonon matrix element , gE 2g = 
Ie 2b (h 2 /4MCI2) 1 / 4 , determines the change in the electron 
energy bands associated with the ground state zero point 
motion. 

Frozen phonon calculations have been previously em- 
ployed in Ref. to point out the connection between 
small Fermi energy and anharmonic effects of the Ei g 
phonon mode in MgB 2 . The appearance of anharmonic 
terms in the static phonon potential has indeed been re- 
lated to a strong lattice displacement regime, namely: 
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FIG. 1: (a) Frozen phonon potential V(u) for the Eig mode in 
MgB2 (solid line, left side scale), and probability distribution 
function P(u) = \^E 2g {u)\ 2 (dashed line, right side scale), u is 
in units of A, V(u) in units of eV and P{u) in units of A~\ (b) 
Electronic structure of undistorted MgB2 (solid lines). The 
grey regions represent the quantum fluctuations induced by 
zero point lattice motion (|tt| < (u 2 ) 1 / 2 ). Bold lines represent 
the two a bands, which are also the most affected by quantum 
fluctuations. 



D E2s = Ie 2!) u > Ep. On physical grounds, we can expect 
that anharmonic terms will be experimentally observable 
only if the lattice fluctuations allowed in the system are 
sufficiently strong to sample regions of phase space in 



which the effective Fermi energy is small (Ie 2 



We are going to show that a similar condition rules the 
breakdown of the adiabatic BO principle. 14 

In Fig. ^a) we plot the E^ g frozen phonon potential 
V(u) obtained by local density (LDA) calculations (see 
Ref. for details). In the spirit of the BO principle 
we can evaluate the ground state phonon wave func- 
tion XE 2g (u) by the numerical solution of the Schrodinger 

equation \-{K 2 /2M E2g )V 2 u + V{u)] X E 2g (u) = E E2gX E 2g . 
We can also define a lattice probability distribution 
functions (PDF) P(u) through the relations P(u) = 
\XE 2g {u)\ 2 , which is plotted in Fig. ^a). In the classical 
(strictly adiabatic) Me 2 — ► 00 limit, and P(u) becomes 
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a (5-function centered at u — 0. The electronic structure 
of the undistorted lattice, which corresponds thus to the 
adiabatic limit, is shown in panel (b) as solid lines. 

As mentioned above, there are two different sources of 
quantum fluctuations of the lattice configuration with re- 
spect to the classical, static limit. The first one is given 
by the possibility of quantum transitions between differ- 
ent electron and lattice eigenfunctions. This gives rise, 
when expanded around the lattice equilibrium position 
u = 0, to the conventional electron-phonon scattering 
which leads to the dynamical electron-phonon renormal- 
ization of the electronic and lattice properties, usually 
taken into account by the corresponding self-energies. A 
different effect of the quantum lattice fluctuations is the 
fact that the ground state PDF acquires a non-zero am- 
plitude at u 7^ 0, corresponding to the zero-point motion. 
In this framework an appropriate quantity to estimate 
the amplitude of the lattice quantum fluctuations is the 
root mean square (rms) lattice displacement: 



.avi/a = 



du u 2 P(u) 



1/2 



(4) 



In the following we focus on the interesting implications 
of this second effect on the electronic structure. We shall 
show that the zero-point lattice fluctuations are unavoid- 
ably reflected, through the deformation potential, in an 
intrinsic source of quantum fluctuations for the electronic 
structure and for the Fermi energy. 

A rough estimate of these effects is provided by their 
rms values. Applying the definition Ijljl. we obtain for 
the zero-point motion of the Ei g phonon mode a vi- 
brational amplitude (it 2 ) 1 / 2 ~ 0.034 A. The changes 
in the electronic structure corresponding to this lattice 
fluctuations are represented in Fig. ^b) by the grey 
regions. The main changes of the electronic structure 
close to the T point, D% 2 — lE 2g {u 2 ) 1 ^ 2 can be ob- 
tained from the deformation potential and roughly corre- 
spond to the thickness of the grey region at the T point. 
As shown in the figure, the strong quantum lattice fluc- 
tuations are reflected in strong electronic changes with 
a rms value of the a band Fermi energy fluctuations 



[Eg (u) - Eg(0)] 
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i 2 ) 1 / 2 ~ 0.39 eV, which 



is slightly smaller than the electron-phonon matrix ele- 
ment gE 2g = 0.45 eV because it includes anharmonic ef- 
fects. This value is comparable to the Fermi energy of the 
undistorted (adiabatic) case Ep = 0.45 eV, resulting in a 
range of Fermi energy fluctuation 0.06 eV < Ep < 0.84 
eV, and casts doubts on the usual definition of the Fermi 
energy, based on the BO approximation. 

The relevance of this regime may be measured by the 
dimensionless parameter: 

: (iEg( M ) -mm 2 ) 1 ' 2 „ i E23 (f) 1/2 ^ 091 (5) 
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In common metals, although the numerator can be of the 
order of a fraction of eV, the denominator is usually of 



5 — 10 eV, so that k <C I and the role of the quantum 
lattice fluctuations on the electronic structure is negligi- 
ble. Things are radically different for the E 2g phonon of 
MgB2, where the large value of k (~ 1) is driven by the 
extremely small Fermi energy of the a bands. 

The sizable magnitude of the parameter k calls for 
some physical considerations. Due to these strong quan- 
tum fluctuations, the system indeed samples, with a siz- 
able weight, electronic configurations with zero or van- 
ishing Fermi energy. The system has thus an intrinsic 
nonadiabatic character— even if the phonon frequencies 
Wph are sensibly smaller than Ep(u = 0). An alterna- 
tive way to understand the breakdown of the BO prin- 
ciple is to observe that for u ~ Ep/lE 2 , the electronic 
Fermi energy Ep(u), which governs the dynamics of the 
a bands, is comparable with w p h, which governs the lat- 
tice dynamics. The variable u cannot thus be treated 
as a classical degree of freedom in a parametric way and 
#(r;u) ^ x(u)(^(r; [u]). 

We would like to stress once more that the origin of the 
nonadiabatic breakdown of the BO principle stems from 
the strong lattice quantum fluctuations which make the 
electron-phonon matrix element, the second term in Eqs. 
@-©i comparable with the Fermi energy. This is thus 
qualitatively different from what discussed in Refs. iTlSl 
where the nonadiabatic effects are ruled by the ratio be- 
tween the energy scale of the third term in Eq. J3J and 
the Fermi energy, namely oje 2 /Ep. Note also that the 
phonon frequency cuE 2g does not depend on the ampli- 
tude of the lattice fluctuations, on the contrary of the 
parameter Ie 29 (u 2 ) 1 ^ 2 /Ep. This means that the onset 
of one kind of nonadiabatic effects does not necessarily 
implies the onset of the other one, and vice versa. 

In order to further clarify this statement let us assume 
for the moment a perfectly harmonic phonon mode with 
clastic constant a 2 . V(u) — a-iv? . Anharmonic contri- 
butions do not play any important role in the following 
discussion, so we will take them into account through 
an "effective" elastic constant as in Ref. [Tj to give the 
anharmonic hardening of the phonon. 

The two different sources of breakdown of the BO prin- 
ciple are illustrated in Fig. [21 where we plot on the x 
axis the inverse of the square root of the phonon mass 
1 / \Z^E 2g and on the y axis the inverse of the square 
root of the elastic constant l/^/a^. The limit Me 2s — >■ oo, 
denoted by the thick line, is the strictly classical (adia- 
batic) case where both kinds of nonadiabatic effects are 
negligible. In the so-defined phase space a fixed phonon 
frequency is given by a straight line cutting through the 
origin. In Fig. [21 we show as dashed lines the three fre- 
quencies u)e 2 = 68 meV, 81 meV, and 0.45 eV, corre- 
sponding to the harmonic and anharmonic frequencies of 
the E2 g phonon, and to a hypothetical phonon in the fully 
nonadiabatic regime (w p h = Ep) respectively. The empty 
(filled) circles represent respectively the actual position 
of the E 2g phonon of MgB 2 (M = 10.81 a.m.u.) with 
(without) taking into account the anharmonic harden- 
ing. The nonadiabatic breakdown of the BO principle 
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FIG. 2: Schematic phase diagram for the different kinds of 
nonadiabatic effects. See text for details. 

underlined in Ref. [7] arises in the region o> p h ~ E^, and 
is represented by the vertically-dashed area, correspond- 
ing to ujph/Ep > 0.5. 

The breakdown of the BO principle due to the zero 
point motion, determined by the parameter re, of the lat- 
tice fluctuation is also shown in Fig. [3 re-isolines are 
easily defined by the relations (^(m 2 ) 1 / 2 = and 
k = Ie 2s (u 2 ) 1/ ' 2 /E f using LDA values I E2g = 12 eV 
/ A and Ep = 0.45 eV. The horizontally-dashed region 
(re > 0.5) represents the range of parameters where nona- 
diabatic effects could be triggered by the zero-point fluc- 
tuations. It is interesting to note that in this new regime 
nonadiabatic effects induced by lattice fluctuations can 
be operative even quite far from the usual "nonadiabatic" 
regime w p h ~ E-p. Of course, the two sources of nonadi- 
abaticity merge together in the highly nonadiabatic case 
M — > 0, where E-p is smaller than w p h already in the 
absence of fluctuations. This could be the case for the 
cuprates and fullcrcnc compounds since in their case Ep 
is smaller than in MgB2, while the relevant phonon fre- 



quencies are comparable. 

In conclusion, in this work we have pointed out the 
inadequacy in MgB2 of the BO approximation which is 
at the basis of many ab-initio techniques, when applied 
to coupling between the Ei g phonon and the a holes, 
which are mainly responsible for superconductivity. We 
have related such inadequacy to the strong lattice fluc- 
tuations of zero-point motion which, due to the large 
deformation potential associated with the E^ g phonon, 
induce electronic changes of the a bands of the same or- 
der of magnitude of their Fermi energy (Ep ~ 0.45 eV). 
The breakdown of the BO approximation has important 
consequences since it implies the failure of many ab-initio 
techniques. An appropriate inclusion of these effects is 
a formidable task which goes beyond the aim of this pa- 
per. We note however that the structural and cohesive 
properties of MgB2 are expected to be weakly affected 
by the zero-point fluctuations, both because they depend 
on the whole set of occupied bands (whose width is 10-20 
eV), and because they involve all phonons. Its transport 
and superconducting properties, instead, are mainly de- 
termined from the small fraction of carriers within an 
energy range ~ w p h close to the Fermi level, and which 
are strongly coupled to the relevant Ei g phonon. 

As a last remark, it is interesting to remind that 
the nonadiabatic effects induced by lattice fluctuations 
are increased as the force constant decreases. Electron- 
phonon systems with incipient lattice instabilities (corre- 
sponding to flattening of the lattice potential) could be 
thus good candidates for the observation of similar nona- 
diabatic effects in other materials than MgB2 . Good can- 
didates in this direction are bismuthates, A15 compounds 
and possibly monoatomic metals as Nb or V. In particu- 
lar, the important role of the lattice fluctuations in A15 
compounds, which was already pointed out by P.W. An- 
derson in Ref. 0, could provide the basis to investigate 
some interesting similarities between A15 and diborides 
systems. 
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